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ABSTRACT 

We report on a sensitive search for redshifted Ha line-emission from three high-metallicity 
damped Lya absorbers (DLAs) at z ~ 2.4 with the Near-infrared Integral Field Spectrometer 
(NIFS) on the Gemini-North telescope, assisted by the ALTtitude conjugate Adaptive optics 
for the InfraRed (ALTAIR) system with a laser guide star. Within the NIFS field-of-view, 
3'.'22 x 2 '!92 corresponding to ps 25 kpc x23 kpc at z = 2.4, we detect no statistically sig¬ 
nificant line-emission at the expected redshifted Ha wavelengths. The measured root-mean- 
square noise fluctuations in 0"4 apertures are 1 — 3 x 10 -18 erg s -1 cm -2 . Our analysis of 
simulated, compact, line-emitting sources yields stringent limits on the star-formation rates 
(SFRs) of the three DLAs, < 2.2 M 0 yr” 1 (3cr) for two absorbers, and < 11 M 0 yr” 1 (3cr) 
for the third, at all impact parameters within pa 12.5 kpc to the quasar sightline at the DLA 
redshift. For the third absorber, the SFR limit is < 4.4 M 0 yr” 1 for locations away from the 
quasar sightline. These results demonstrate the potential of adaptive optics-assisted, integral 
field unit searches for galaxies associated with high -z DLAs. 

Key words: quasars: absorption lines - galaxies: high redshift - quasars: individual: 
Q0201+365 - quasars: individual: Q0311+430 - quasars: individual: Q2343+125 


1 INTRODUCTION 

The damped Lya systems (DLAs: Iwolfe et al J1 19861 . 120051) de¬ 
fine the class of absorption-line systems discovered in the rest- 
frame UV spectra of distant quasars, with HI column densities 
JVhi ^ 2 x 10 20 cm” 2 as measured from the analysis of damping 
wings in the Lya profile. These large /Vhi values, characteristic 
of gas in th e interstellar m edium of modern, star-forming galax¬ 
ies (e.g. lWalter et ahlfcoosl) . have led researchers to associate high- 

2 DLAs with young galaxies in the early Universe. This associa¬ 
tion is well supported by theoretical models of galaxy form ation 
(e.g. lHaehnelt et al.ll 19981 : [Mailer et alJl200ll:lBird et alj|2014h and 
recent numerical simulations that predict large HI surface densi¬ 
ties extending out to tens of kpc from active, sta r-forming regions 
(e.g. lRahmati et al.ll2013l : lFumagalli et al.ll2014al) and even to dis- 
tances of pa 200 kpc from the centres of the most massive halos 
JPontzen et alfcoda ). Further empirical evidence connecting high- 
2 DLAs to galaxies includes the omnipresence of heavy elements in 
the gas typically at levels exceeding estimates for the intergalactic 
medium which indicates previous, if not recent, pollution by stars 
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(e.g.lPetti ni et al.ll 994lll997l : |Prochaska et alj2003l : |Penprase et al.1 


l2010l:lRafelski et alj|2012k At low redshifts, where one may more 

sensitively search for associated galaxies, sever al st u dies have id en¬ 
tified galaxy counterparts to the a bsorbers (e.g. lle Brun et all 19971 : 
iRao et aljEzOollChen et all2005t) . 

Detailed analysis of DLA absorption lines offer unparalleled 
insight into the interstellar medium (ISM) of these galaxies, 
including their chemical abundance patterns, gas temperature. 


molec ula r content, dust depleti on, and k ine matics (e.g. Pettini et al 


2003 


1997; Prochaska & Wolfel 19971: Kanekar & Chengalui 


|Ledoux et al. |2003|: |Dessauges- Zavadskv et a lJ 120061 : 


Noterdaeme et al 20081 : 1 Kanekar et alj2014l) . Indeed. DLA studies 
provided the first and (still) the most comprehensive view of 
the ISM at high redshifts. To properly place this wealth of ISM 
measurements within galaxy formation theories, however, it is 
critical to connect the gas to the stars and to the processes of star 
formation. Establishing this connection directly by observing the 
stars and/or HII regions, however, has proven to be remarkably 
challenging. 

Early works searched first for Lya emission, both in 
the core of the HI Lya absorption profile and with narrow- 
band imaging_or, more recently, with integra l field spectroscopy 
(e.g. ISmith et al.l 1 19891: iHunstead et aD Il990l : IWolfe etalj 1 1992: 
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M 0 II 1 


er & Warrenl 1 19931 : iKulkami et al.l 1 20061 : |Christensen_rtal 


2007). These provided a few rare detections le.g. lDiorgovski et al 


19961 : iFvnbo et alii 9991 : iMelller et al.ll2004l) and mostly upper lim¬ 


its to the star-formation rates (SFRs), 10 M@ yr _ , albeit sub¬ 
ject to substantial systematic uncertainties related to dust depletion 
and Lya radiative transfer. The community then searched sensi¬ 
tively for far-UV continuum emission in broad-band imaging, ap¬ 
plying the Lyman-break tech nique when possible to pre-select can¬ 
didate s at high redshifts (e.g. iMoller & Warrenll 19931 : IWarren et al .1 
l200ll : IProchaska et al J 12002) : these surveys were largely unsuc¬ 
cessful, with only a few positive detections. This implied inher¬ 
ently faint ultraviolet emission and/or that the absorbing galaxies 
were drowned out by the glare of the background quasar. More 
recently, IFumagalli et al J ( 120100 introduced a novel technique to 
address the problem of detecting a faint DLA close to a bright 
quasar, using a high-iVHi absorption system at a redshift close to 
that of the background quasar to “block” its far-ultraviolet light, 
thus allowing a sensitive search for emission from a second, lower- 
redshift DLA along the sightline. Unfortunately, even this tech¬ 
nique has so far yielded null results, albeit yielding more stringent 
constraints, SFR ^ 0.1 — 0.3 Mq yr _1 at all impact parameters 
dFumagalli et all2014bl , [20l3) . 

Somewhat stymied, the community has made two changes in 
the observational strategy: (1) focus the effort on the most heav¬ 
ily enriched DLAs under the expectation that these will be as¬ 
sociated with bright er galaxies having higher SFRs and/or larger 
stellar masses (e.g. iMbller et~al1 |2004 ILedoux et al.l 1 20061) . and 
(2) search for the strong nebular emission lines associated with H II 
regions (e.g. OIII, Ha). Earlier studies using the latter approach 
yielded relatively weak upper limits on the SFRs of DLAs, < 10 — 
100 M w yr -1 at z }' 2 (e.g.lBec htold et a lJI 19981 : IMannucci et al.l 
1 19981: iBunker et al.ll 19991: iKulkarni et alj|200d~ when converted to 
the LCDM cosmology used in this paper). However, this method 
has recently been significantly advanced by technical progress, i.e. 
the commissioning of higher throughput near-infrared (near-IR) 
spectrometers on large-aperture telescopes. Several of these include 
integral field units (IFUs) which permit a “blind” search for emis¬ 
sion in a modest field around each DLA. These new efforts, mostly 
with the Very Large Telescope (VLT), have met with somewhat 
greater success: ~ 10 counterparts to DLAs at z > 2 have recently 
been discovered, at impact parameters from the quasar sightlines 
in th e range ~ 1 — 23kptQ and with SFRs 10 — 20 Mq yr~ 


(e.g. Fvnbo e£al 


Noterdaeme et al 


20 l’(j 1201 lL 120131: P eroux et al. 1201 ll. l20ll 
120121) . There have also been a few recent detec¬ 


tions ofL^o_et^smnin thetrough_oftheDLAL^oabsorgtionJhne 


(e.g.lMoller et al.ll2004l:lHennawi et alEoollKulkami et alj|2012l: 


iKrogager et all2012l:lNoterdaeme et al.ll2014l). 

One is scientifically motivated to search for star-formation 
from DLAs for several reasons beyond mere studies of this ab¬ 
sorber population. In contrast t o the multitude of surveys of star¬ 
forming galaxies at z ~ 2 (e .g. iDaddi et aPl2004l : lErb et al.ll2006l : 
iForster Schreiber et al.ll2009l) . one may link the gas directly to the 
stars and star-forming regions in these young galaxies. For HI 
in particular, absorption-line analysis represents the only means 
for assessing the surface density of gas around these galaxies and 
thereby an estimate of the gaseous extent. 

Studies of the star-forming regions of DLAs also offer a com- 


1 Throughout this paper, we will use a flat A-cold dark matter cosmol¬ 
ogy, with Qm = 0.315, Da = 0.685 and Ho = 67.3 km s~ 1 Mpc —1 
i Planck Collaborati o n |201~3) . 


plementary assessment of emerging trends in z ~ 2 galaxies be¬ 
tween stellar mass M*, metallicity and star-formation rate. Large 
spectroscopic surveys have revealed that star-forming galaxies fol¬ 
low a relatively tight main sequence, in which the SFR correlates 
with the stellar mass ( e.g. iNoeske et al]|2007l : IDaddi et all 1 20071 : 
iRodighiero et al][201 lh . This has inspired a series of simple mod¬ 
els of galaxy formation within the ACDM framework that con¬ 
sider the balance of gas accretion, star formation and the_out- 
flow of materia l by feedback pro cesses (e.g. iDekel et al] 120091 : 
Lilly et al1 i2013l : lForbes et al.l2014l) . Observers have also explored 
second-order relations in the galaxy properties, emphasizing the 
existence of a fundamental metallicity relation (FMR) where the 
SFR is sys tematically higher in lower metallicity galaxies at a 
fixed mass ( IMannucci et al.ll20ld : lBothwell et aT]|20 1 3l : [Stott et all 
120131) . With the DLAs, one may examine whether absorption- 
selected galaxies follow the same trends as these traditional sam¬ 
ples. In addition, one may further test the simple models by incor¬ 
porating the ISM properties constrained directly by DLA observa¬ 
tions (e.g. HI extent, ISM metallicity, star formation history, gas 
temperature, etc). 

In this manuscript, we report on our own project to search for 
Ha emission from high-z DLAs_with the Near-infrared Integral 
Field Spectrometer (NIFS; iMcGregor et al.ll2003h on the Gernini- 
North telescope. We pre-selected a set of DLAs with high metal¬ 
licity and with absorption redshifts well-suited to the detection of 
nebular emission lines (e.g. strong sky lines were avoided). In ad¬ 
dition, we used the ALTtitude conjugate Adaptive optics for the 
InfraRed (ALTAIR) with a laser guide star to isolate the light of the 
background quasar and maximize the signal from these (presum- 
ably) small galaxies. We note red that I Jorgenson & Wolfel d2014l) 
also employed a similar approach, with laser-guided adaptive op¬ 
tics and an IFU with the Keck I telescope to analyze a previously 
known galaxy counterpart toaz ~ 2.3543 DLA with [O III] and 
Ha emission. 

This paper is organized as follows: §[2]describes the DLA sam¬ 
ple, the observational programme, and the details of the data reduc¬ 
tion. The primary results are given in §[3]and we conclude with a 
discussion of our non-detections, similar searches in the literature, 
and future efforts in § [4] 


2 OBSERVATIONS AND DATA REDUCTION 

From the set of known DL As at z ~ 2 in 2007 , with accurately 
meas ured metallicities (e.g. ILedoux et al J1 20061 : IProchaska et al.l 
120071) , we considered all systems with enrichment exceeding ~ 
1/3 solar abundance. We further restricted this set to have good 
visibility from Mauna Kea and chose five systems with absorption 
redshifts z a bs that placed their Ha lines at favorable wavelengths 
for deep, near-IR spectroscopy (i.e. avoiding bright sky lines and 
strong telluric absorption). Unfortunately, two of the DLAs were 
observed with incorrect wavelength settings, so we were left with 
three absorbers; TableQ]summarizes their properties. 

We used NIFS on Gemini-North to search for redshifted Ha 
emission from the three DLAs in 2008, with ALTAIR and laser 
guide stars to improve the point spread function (PSF). Nearby stars 
or the QSO itself were used for the tip-tilt mode correction in the 
laser guiding. The data were acquired in queue-observing mode 
during semester 2008B under Program ID GN-2008B-Q-61. The 
observing conditions were photometric with natural seeing better 
than 0"65 in the optical V-band as measured by the ALTAIR wave- 
front sensor. 
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Ha emission in high-z DLAs 3 


Table 1 . DLA Sample 


QSO 

RA (J2000) 

DEC (J2000) 

^em 

2DLA 

JVhi 

[Z/Hp 

References 

QSO0201+365 

02:04:55.60 

+36:49:18 

2.912 

2.4628 

20.38 

-0.29 

1,2 

QSO0311+430 

03:14:43.60 

+43:14:05.1 

2.870 

2.2898 

20.30 

-0.49 

3,4,5 

QS02343+125 

23:46:28.22 

+12:48:59.9 

2.763 

2.4313 

20.34 

-0.54+ 

6, 7,8 


Notes to the table: * The metallicity estimates are based on [Zn/H] (QS00201+365) and [Si/H] (QSO0311+430 and QS02343+125). 
t The metallicity of the z = 2.4313 DLA towards QS02343+125 is from the HIR ES spectrum o f Lu et all J19981) . We have re-analysed these data and obtain 


a slightly higher, albeit consi stent, metallicity of [Si/H]= —0.4 7. We note that lLedoux et alj 


20061) obtained a lower metallicity from the Znll lines. 


[Zn/H]= —0.92 ± 0.07, while iDessauges-Zavadskv et alj f2004h u sed Sll lines to o btain [S/H]= —0.7 ±0.1 (consistent with the metallicity estimate of 

___ i _ iLu et akU 19981)). ____ ___ ( ( _ 

References: (l) IProchaska & Wolfel 19961 : (2) |Procha ska et al l 200 it (3) York et al j2007[: (4)lEllison et al.l200d (5) lKanekar et all2Q14l : (6) lLu et all 19981 : 

(7)lDessauges-Zavadsky et a ll2004 f8ilProchaska et al.ll2007t 


We used the K grating to search for Ha emission from the 
three DLAs at z « 2.4. For each target, the observing band cov¬ 
ered a wavelength range of approximately 4000 A, centered at the 
expected DLA Ha wavelength (see Table [2j. The NIFS field-of- 
view is approximately 3” on a side, corresponding to w 25 kpc at 
z = 2.4. Each exposure was set to 560 seconds, a compromise 
between minimizing the effects of instrumental flexure and read¬ 
out noise. The exposures were dithered along R.A. and Dec. with 
±0715 offsets to obtain better estimates of the sky background. 
This offset was chosen to match the expected PSF of ~ 07 15. How¬ 
ever, the actually achieved PSF is worse than this (see below), lead¬ 
ing to a less optimal sky subtraction for slightly extended objects. 
We did not take off-target exposures to independently measure the 
sky background. A0V stars at similar airmass were observed imme¬ 
diately before and/or after the target observations with an identical 
instrumental setup to estimate telluric absorption and to perform 
flux calibration. Images of Ar and Xe arc lamps were obtained to 
provide wavelength calibration. 

The data were reduced with our own programmes written in 
the Interactive Data Language environment. The images were first 
dark subtracted and flattened; Figure |T}i shows an example of one 
processed frame. As is evident from the figure, the NIFS lenslets 
project spectra onto the detector with different wavelength zero- 
points. To begin sky subtraction, we first registered the sky lines 
by applying a rigid pixel offset to each row derived from a 1-D 
cross-correlation, after a 2 x re-sampling in the spectral dimension. 
Because the 2x re-sampled spectral pixels have a size of 1.03 A, 
this rigid pixel offset does not introduce a wavelength error larger 
than 0.5 A. 

After the rows of spectra were aligned, a constant sky back¬ 
ground was subtracted column by column (Figure|T[j). At this stage, 
this simple sky-subtraction was sufficient to allow detection of the 
QSOs over the entire spectral range, in each 560s exposure. The 
original NIFS 2-D images have a spatial pixel size of 07103 x 0704 
and spectral pixel sizes of 2.07 A. The images were reformatted and 
re-sampled into 3-D cubes of 0 , /052 x 0704 x 2.07 A pixel sizes. 
The cubes were collapsed along the spectral axes to form contin¬ 
uum images of the QSOs (Figure [2] top panels). Centroids of the 
quasar emission were used to derive dither offsets for the multiple 
exposures with a precision that is much smaller than the pixel size. 
We also estimated the PSF provided by the ALTAIR adaptive op¬ 
tics (AO) system to lie in the range 0.2 — 0.3" (full width at half 
maximum; FWHM). 

To sensitively search for the DLA Ha emission, one requires a 
much higher fidelity approach to sky subtraction and cosmic ray re¬ 
moval. FigureQJ) shows obvious residuals at the brightest sky lines 
and also “ripples” in the spatial dimension that we associate with 



Figure 1 . Example of our NIFS data processing (see main text for details). 
In all panels, the horizontal and vertical directions are the spectral and spa¬ 
tial directions, respectively, (a) shows an image after dark subtraction and 
flat fielding, (b) shows the result after the vertical sky lines are horizontally 
aligned and subtracted. The bright horizontal lines near the middle of the 
image are the QSO continuum, and the broader horizontal bands are resid¬ 
ual sky background, (c) shows the result after the subtraction of a smooth 
function, column by column and then row by row. (d) shows the result after 
cosmic ray removal and subtraction of an image that is median-combined 
from all dithered exposures. 


spatial variations in the sky background and/or non-uniformities in 
the lenslets that have not been captured by the flat-fielding proce¬ 
dure. We also wished to search sensitively for Ha emission close 
to, or even at, the position of each QSO. Therefore, we fitted a 
5th or 6th-order polynomial, column by column and then row-by- 
row, to each NIFS 2D image and subtracted the fitted values. Fig¬ 
ure m shows a resultant image. This procedure is designed to re¬ 
move slowly varying flux fluctuations related to the background 
and the QSO without eliminating narrow features such as emis¬ 
sion lines and cosmic rays. For the latter, we flagged all pixels with 
fluxes 8a above the background and verified (by eye) that none 
of these could be consistent with an emission feature. Lastly, we 
median-combined the five to ten dithered images without alignment 
to generate a final residual map of the sky background including 
the effects of “hot” pixels. This residual map was subtracted from 
each image (e.g. Figure |TJl) and these images were reformatted to 
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Table 2. Observing Logs 


Field 

Date (UT) 

Kbs (pm) 

Telluric Stars 

^"exp (sec) 

Magtt 

dtt 

QS02343+125 

2008 Oct 4 

2.2519 

HD 21501 

4480 

V = 17.0 

0" 

QS00201+365 

2008 Oct 7 

2.2726 

HD 1561, HD 21038 

5600 

R = 16.2 

1970 

QSO0311+430 

2008 Oct 7 

2.1590 

HD 12559 

2800 

F = 16.0 

2674 


Magtt and dtt are the magnitudes of the tip-tilt guide stars and their angular distances to the targets, respectively. 




QSO Continuum 


QSO/Sky Subtracted 
300 km/s Average 


QSO/Sky Subtracted 
600 km/s Average 




Figure 2. Results of our Gemini-NIFS search for Ha emission from three 
metal-rich DLAs. Each panel has an angular size of 3722 X 2792 (corre¬ 
sponding to ~ 25kpc X 23 kpc at z = 2.4) and has an inverted gray scale. 
The top row shows continuum images of the three background QSOs. The 
middle and bottom rows show averaged images within 10 and 22 spectral 
pixels (corresponding to ~ 300 km s~ 1 and ~ 600 km s” 1 ), respectively, 
centered at the expected redshifted Ha wavelengths, after subtraction of 
the QSOs and the sky background. No evidence of Ha emission is found in 
these images or in the data cubes. 


3-D cubes. The dithered cubes of each target were then aligned and 
combined with sigma-clipping to form the final data cube for each 
field. 

The spectral images of the Ar/Xe arc lamps and the A0V 
stars were processed separately. Sky-subtraction was not carried 
out here since the sources are much brighter than the structure in 
the sky background. The extracted Ar/Xe spectra were used to de¬ 
rive wavelength solutions, accurate to ~ 1 A over the observed 
spectral ranges, which is sufficient for our science. The extracted 
A0V spectra were used to derive telluric absorptions and perform 
flux calibration. First, we masked the 2.166 gm hydrogen absorp¬ 
tion lines on the A0V spectra and interpolated over the masked 
regions in the continua. T he spectra w ere then compared with the 
model spectrum of Vega iKuruczlI 19791) to derive the telluric ab¬ 
sorptions as a function of wavelength. Absorption corrections of 
« 1% were applied to both the data cubes of the targets and the 
A0V spectra. Finally, the absolute flux calibrations were estimated 
with the 2MASS Ks fluxes of the A0V stars and their observed 
spectral fluxes in the 2MASS Ks passband (~ 2.17 /im). The flux 
zero points derived from the four A0V stars under three different 
spectral settings agree very well with each other (to within ~ 3%), 
suggesting that our flux calibration is accurate to within 5%. 




Figure 3. Histograms of the fluxes measured in 072 (dotted), 074 (solid), 
and 076 (dashed) square apertures across the images in each field. The mea¬ 
surements were made within a ~ 300 km s —1 velocity range, centered on 
the expected redshifted Ha wavelength. 


3 RESULTS 

The final reduced images are shown in Figure [2] Each 
panel is 3722 x 2792 in size, corresponding to approximately 
25 kpc x 23 kpc at z = 2.4. The top row shows the stacked 
QSO continuum images from an intermediate stage before pro¬ 
cessing (see Here, the sky subtraction is very simple and 
residual sky structures are apparent. In the case of QS00201+365 
and QSO0311+430, the tip-tilt guide stars are quite far away 
from the QSO and we only obtained partially corrected PSFs 
from ALTAIR. The measured PSF FWHMs are O'.'24 x 0'.'17 for 
QS02343+125, 0736 x 0723 for QS00201+365, and O'.'31 x O'.'23 
for QSO0311+430. These are all much better than the natural see¬ 
ing but worse than what ALTAIR can achieve in the best cases 
(~ 07l). The asymmetric PSFs are likely to arise from the laser- 
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guided AO performance, as the self-guided observations of the 
AOV stars do not show such asymmetry. 

The middle row of the figure shows stacked images spanning 
10 spectral pixels (corresponding to ~ 300 km s -1 ) centered at 
the expected redshifted Ha wavelength for each DLA, after fur¬ 
ther subtraction of the QSO and sky background (see §0. Both the 
QSO and the sky were cleanly subtracted across the images, en¬ 
abling the search for faint Ha emission even along the sightline to 
each QSO. There is no evidence of significant Ha emission in any 
of the three DLAs. We also carefully searched the 3-D image cubes 
over a broader wavelength range, in images before and after cos¬ 
mic ray removal, and in individual unstacked exposures. We found 
no evidence for emission in any of these cases. To search for very 
broad emission, we also created images averaged over 22 spectral 
pixels (corresponding to ~ 600 km s -1 , bottom row of Figure [2jt; 
no emission was detected at this coarser resolution. For the discus¬ 
sion below, we will assume a line width of ~ 300 km s , i.e. we 
will use the images of the middle row in Fig. [2] 

During the data reduction process, we median-combined the 
dithered exposures to form residual images for better sky subtrac¬ 
tion, using dither offsets between exposures of « 0” 15. This means 
that the DLA Ha emission might be partially subtracted if it is 
smooth and extended over angular scales of > '.'3. We did not 
see any hint of such extended emission in the images before the 
median sky subtraction, nor in the median-combined sky image. 
Furthermore, such spatially-extended, but spectrally-narrow, emis¬ 
sion ( exten ded over > 2.5 kpc at z = 2.4) appears unlikely (e.g. 
iFvnbo et al.l2oTHl . We conclude that the three DLAs of our sample 
do not produce strong Ha emission within ~ 12.5 kpc of the QSO 
sightline. 

We used the non-detections of Ha emission to estimate up¬ 
per limits to the Ha flux from each DLA. We measured the back¬ 
ground flux fluctuations in the 300 km s -1 images shown in the 
middle row of Figure [2] in (//4 square apertures across each im¬ 
age. The O'/4 size is larger than the PSF FWHM in all three im¬ 
ages, and corresponds to 3.3 kpc at z = 2.4. Figure [3] shows 
a histogram of these fluxes, as well as results for other aper¬ 
ture sizes. The measured root-mean-square (RMS) noise fluctua- 

1.38 x 


tions in O'/4 apertures are 2.86 x 10 18 erg s 1 cm 2 
10 - 


velocity profiles over 10 —15 A. (« 100 — 200 km sec -1 ) and were 
randomly placed in the raw images with a range of fluxes. We con¬ 
sider a simulated source as recovered if it passes our inspection and 
its flux is detected at > 3a significance. We found that we con¬ 
fidently recover the simulated sources in our reduction when the 
input flux is ~ 4 times larger than the formal RMS noise. This 
is consistent with the above-mentioned 20%-30% flux loss caused 
by the broad PSF. We thus regard this as a 3a confidence detec¬ 
tion limit. Note that this is true even when the simulated sources 
have very small impact parameters to the QSO sightline, except 
for the case of QS02343+125. QS02343+125 has broad emission 
lines in our observed spectral cube, making our quasar subtrac¬ 
tion slightly poorer. In this case, we can only detect the simulated 
line source along the QSO sightline when it has > 10a signifi¬ 
cance. We therefore conclude that our observations constrain the 
star formation rates to SFR < 2.2 M@ yr -1 at all impact parame¬ 
ters for the DLAs toward QS00201+365 and QSO0311+430, and 
SFR < 4.4 M@ yr -1 at separations > O'/5 for the DLA towards 
QS02343+125, with a poorer limit of SFR < 11 M 0 yr -1 close 
to the QSO for the last system. We note that all of these limits ig¬ 
nore dust obscuration be yond th e m odest correction associated with 
the empirical relation of lKennicutd Il998l) . Our observations would 
thus underestimate the SFR limits in the case of modest amount of 
extinction, and would be insensitive to dusty star-forming regions, 
for which the Ha line emission might be highly obscured. 


4 DISCUSSION 

The primary aim of searches for line emission from high-3 DLAs 
is the direct identification of their galaxy counterparts, so as to es¬ 
timate the typical size, mass and SFR of absorption-selected galax¬ 
ies. Non-detections of Ha emission only allow us to place limits 
on the SFR, assuming that the emission region of the DLA host 
lies within the search area. In addition, however, there has been 
much recent interest in whether high-3 absorption-selected galax¬ 
ies also have a mass-metallicity rel ation similar to that se e n in high- 
z em ission-selected galaxies (e.g. iTremonti et all 1 20041 : lErb et al,l 
1 20061) . The existence of such a relation in DLAs is suggested 


in the 


images for QS02343+125, 

QS00201+365, and QSO0311+430, 

velocity width of low-ionization metal lines ( 

Wolfe & Prochaska 

respectively (see Figure |3V 

The above values can be converted 

1998; Ledoux et alj2006i;lNeeleman et al. 201; 

) and the rest equiv- 

to limits on the SFR via 

the expression SFR = 7.9 x 

alent width of the Si IlA 1526 line (Prochaska et al. 200811. Simi- 


(erg s -1 ) -1 


Cenmcutt 


10 I/Ha M 0 yr - 
Salpeter initial mass function is assumed. The results are la SFR 
limits of 1.11 M 0 yr -1 , 0.55 M 0 yr -1 , and 0.53 M 0 yr -1 , for 
QS02343+125, QS00201+365, and QSO0311+430, respectively. 
For a line width of 600 km s -1 , the RMS noise and the SFR limits 
increase approximately by a factor of v/2. 

Our relatively small dither offsets (O'/15) and the non-ideal 
PSFs lead to 20%-30% (depending on the source size and PSF) 
flux losses in the median sky subtraction. As a result, the sensi¬ 
tivity to Ha line emission is somewhat worse than the typically 
quoted 3a upper limits based on the above la error estimates. In 
order to better understand the systematics arising from our observ¬ 
ing and data analysis procedure, and to obtain reliable constraints 
on the SFRs, we created artificial line emitting sources in each of 
the raw images of our target sources and processed these images 
identically to the actual images. We assumed intrinsically extended 
emitters with various sizes of Gaussians of up to O'/2 FWHM. We 
then convolved the source images with the PSF derived from the 
QSO images. The sources were conservatively assumed to have flat 


larly, the anti-correlati on detected between DLA metallicity and 
gas spin temperatu re l lKanekar et alj 1 20091 : lEllison et alj 120121 : 
IKanekar et al.l 120141) also supports such a mass-metallicity rela¬ 
tion, as higher mass DLAs are expected to ha ve larger fractions of 
cold g as and hence, lower spin temperatures (IKanekar & Chengalud 


12001 ). To this end, even non-detections of Ha emission in high- 


metallicity DLAs are of much interest, as these may argue against 
the existence of a mass-metallicity relation in D LAs giv en the sta r- 
forming main sequence of 3 ~ 2 galaxies (e.g. lDaddi et al.ll2007lf . 

Before we consider the repercussions of our results, we briefly 
discuss other recent searches for line emission from high-3 DLAs. 
There have been two broad strategies that have been successful in 
either detecting, or placing strong const raints on, Lya, Ha, Oil 
or O III line emission. (1) IFvnbo et alj d2010h pioneered an ap¬ 
proach based on using VLT X-shooter slit spectroscopy, cover¬ 
ing a large wavelength range (3100 A —2.5/rm) to simultaneously 
search for a number of lines, an d using 3 slit position angles centred 


on the QSO location (s ee also IFvnbo et alj 201 ll: Krogager et al.l 


1 20 12l : IFvnbo et al.|[2013lf . (2) Conversely. IPeroux et al.1 d201 il l used 
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IFU spectroscopy with the Spectrograph for Integral Field Ob¬ 
servations in the Near-Infrared (SINFONI) onboard the VLT to 
search for redshifte d Ha emission, simi lar to our approach (see also 
IBouche et al.ll2012l ; IPeroux et alj|2012h . Fynbo et al. also targeted 
high-metallicity DLAs at z > 2, while Peroux et al. and Bouche 
et al. targeted DLAs, sub-DLAs and strong Mg II absorbers of all 
metallicities, centred at two redshifts, 2 ~ 1 and z ~ 2. 

There are pros and cons to the two literature approaches, as 
well as to our ow n approach. The X-shooter-based strategy of 
iFvnbo et al] J2010h has the advantage that a number of lines are 
simultaneously targeted, implying fewer problems due to night sky 
lines and sky subtraction. Also, the use of 3 slit angles reduces the 
likelihood that the DLA might be missed by the slit. However, the 
likelihood of missing the DLA increases with distance from the 
QSO location (as the slits are centred on the QSO), implying a bias 
towards absorbers at low impact parameter, < l". 

The VLT-SINFONI IFU searches for Ha emission by 
IPeroux et alj j 201 2l ) and IBouche et alj j 201 2l) have the advantage 
of being sensitive to Ha emission from a relatively large area (field 
of view ~ 10 x 10 arcsec 2 ) around the QSO sightline. However, 
AO was used (with nearby guide stars) for only a few DLAs and 
even these only yielded a seeing ^ (// 5. The relatively poor seeing 
and the degradation in the signal-to-noise ratio at the QSO location 
(due to the higher background) implies that the sensitivity of such 
searches is lower close to the QSO. Non-detections of Ha emission 
might hence arise a t low impact parameters, < O'/5 (although see 
Bouche et afl llO 1 2l) . 

In the case of our own search, the use of AO implies excel¬ 
lent seeing and high sensitivity to Ha emission even at the QSO 
location for two of the targets, and reasonable sensitivity for the 
third target. However, the relatively small NIFS field of view im¬ 
plies that we are only sensitive to emission for impact parame¬ 
ters < 12.5 kpc. Thus, Ha emission from large disks (of radius 
10 kpc), where the Lya absorption occurs in the outskirts of the 
galaxy while the Ha emission is dominated by the central region 
might be missed in our search. Note that all three DLAs of our sam¬ 
ple have relatively low Hi column densities, < 3 x 10 20 cm 2 , as 
might be expected in the outer regions of a galaxy. 

Moving to detections, and focusing on DLAs at 2 > 2 (i.e. ex¬ 
cluding low -2 DLAs and Mg II absorbers without IVhi estimates), 
the X-shooter method has led to the detection of line e mission 
from five high-meta l licity DLAs at 2 > 2 l lFvnbo et alJl201(iL 
IKrogager et al.ll2012l ; IFvnbo et al.ll20l 3l) . Non-detections from this 
approach have so far not bee n reported. Perhaps u nsurprisingly, 
four of the five detections (see IKrogager et aPl2012h have been at 
relatively low impact parameters, < Cf/S. _ 

In the case of VLT-SINFONI IFU spectroscopy, IPeroux et al .1 
j2012l) report only a single detection of Ha emission in 10 DLAs 
and 2 sub-DLAs at 2 > 2. Further, even the sole DLA of their 
sample with a detection of Ha emission (at 2 = 2.3543 to¬ 
wards J2222—0946) had earlier been detected in Ho, Lya OII and 
OIII lines in X-shooter spectroscopy jFvnbo et alj2Qlof) . This sys¬ 
tem also has an intermediate impact parameter, b « 0 / /8, larger 
than the seeing du ring the VLT-SINFONI observati ons (~ O'.'fi; 
IPeroux et alj|2012h . Note that only two DLAs of the IPeroux et al] 
d2012h sample are “high-metallicity” systems, with [Z/H] > —0.5; 
one of these is the absorber towards J2222—0946, while the other, 
at 2 = 2.2100 towards QSO2059—0528, yielded an Ha non¬ 
detection, with an upper limit to the SFR, < 3M 0 yr 1 0 


2 Note that different authors have used different spatial and spectral resolu- 


iBouche et alj ( l2012h also used VLT-SINFONI to search for 
Ha emission in a set of absorbers at 2 > 2, selected either due 
to strong Mg II absorption or due to identification as DLAs or sub- 
DLAs. We will here only consider the subset of 10 absorbers, 5 
DLAs and 5 sub-DLAs, with Nm estimates. These yielded two 
detections of Ha emission, one system again being the ubiquitous 
2 = 2.3543 DLA towards J2222—0946. The sole other detection 
was in the 2 = 2.3288 DLA towards HE2243—6031, an interme¬ 
diate metallicity absorber ([Zn/H] = —1.10 ± 0.05; iLopez et al .1 
120021) . detected at a large im pact parameter, ~ 2 6.5 kpc, with an 
SFR of 17M 0 yr -1 (see also lBouche et al.ll2013l) . The seven non- 
detections yielded typical SFR limits of 3.2 M 0 yr -1 . 

Finally, our own search in three DLAs yielded no Ha detec¬ 
tions, with strong SFR constraints of 2.2 — T.TMoyr” 1 within 
a distance of 12.5 kpc to the QSO sightline (excepting the QSO 
location in one field). All three DLAs have high metallicities, 
[Z/H] > —0.54, within the top 10% of DLAs at all redshifts (e.g. 
iRafelski et al.ll2012l) . 

It thus appears that, despite a number of studies with 8m- 
class telescopes, there have been very few IFU detections of Ha 
emission in DLAs at 2 > 2. Speci ficall y, combining our re - 
sults with those of IPeroux et alj J2012I) and IBouche et alj J2012I) . 
there has been only a single “new” detection of Ha emission, 
and a corresponding identification of a galaxy counterpart, and 
21 non-detections from IFU searches. Conversely, there have been 
five identifications of g alaxy counterparts via the X-shooter slit 
spect roscopy approach iFvnbo et al] I2010L 1201 ll ; IKrogager et al] 
l2012l) . The target DLA samples of the two types of searches have 
been very different, with the IFU searches mostly targeting low- 
metallicity DLAs (18/22) and the X-shooter spectroscopy focus¬ 
ing on high-metallicity DLAs. This suggests that high-metallicity 
DLAs have higher SFRs, as w ould be expected from a mass- 
metal l icity relation in DL As (e.g. lMetller et alJl2004l ; lLedoux et all 
120061; M 0lleretalJl2Ol3) . together with the observed SFR-M* 
main sequence. 

On the other hand, the three high-metallicity DLAs of our 
sample were not detected in a sensitive IFU search. These have 
metallicities s imilar to those of three of the X-shooter targets in 
IKrogager et alj j2012h . and would have been detected in our search 
for similar SFRs (> 2.2 — 4.4M 0 yr” 1 ) a nd impact parameter s (< 
8 kpc). Adopting the SFR-M* relation of lDaddi et al.1 d2007l) and 
extrapolating it to lower SFRs, our limits imply M* < 1O 9,1 M 0 . 
Is such a low mas s inconsistent with the DLA metallicities? Evalu¬ 
ating the FMR of lMannucci et al.l j20ld) . which is based on emis¬ 
sion lines, with M* = 1O 91 M 0 and SFR = 2.2M 0 yr” 1 , we 
recover 12+log(0/H) « 8.3 dex. This value is roughly consistent 
with the DLA absorption metallicity, especially given the system¬ 
atic uncertainty of nebular-line metallicity measurements inherent 
to the FMR and the possibility of abundance gradients. 

Of course, the number of high-metallicity DLAs that have 
been searched for Ha or nebular emission is still quite small and 
it is possible that observational biases might have caused som e 
of the non-detections. For example, in the model of IFvnbo et al .1 
0008 ). high-metallicity DLAs are expected to arise in large, lu¬ 
minous galaxies, with statistically larger impact parameters. The 
Ha non-detections in our searches might then have been due to the 


tions, as well as significances, when quoting upper limits. Here, we adjusted 
their values to a fiducial 3cr significance and have converted to the same 
Ha calibration adopted here. Differences between the adopted cosmologies 
have negligible impact on the results. 
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limited NIFS field of view. We conclude that, although the paucity 
of detections in IFU searches is certainly curious, all three ap¬ 
proaches (VLT-X-shooter, Gemini-NIFS and VLT-SINFONI) have 
biases (see above) that would affect the fraction of detections of Hct 
or nebular emission, and hence the inferred SFR limits. We suggest 
that the present approaches should be augmented by (1) combin¬ 
ing wide-held IFU searches with the use of AO with laser guide 
stars, so that Ha emission cannot be hidden under the QSO emis¬ 
sion or missed by being outside the held of view, and (2) following 
up X-shooter non-detections with more slit angles, to rule out the 
possibility that non-detections arise because the nebular emission 
is not covered by the slit. Finally, given the biases in the current 
searches, we also feel that it is premature to interpret the present 
SFR data to either support or rule out a mass-metallicity relation in 
damped Lya systems. 
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